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Abstract. We report on the performance of Lagrangian perturbation 
theory up to the second order for the standard cold dark matter (SCDM) 
and broken scale invariance (BSI) scenarios. We normalize both models 
to the COBE data, the BSI model serves as an example of models which 
fit the small-scale power of galaxy surveys. We optimize Lagrangian 
perturbation solutions by removing small-scale power from the initial data 
and compare the results with those of numerical simulations. We find 
an excellent performance of the optimized Lagrangian schemes down to 
scales around the correlation length or smaller, depending on the statistics 
used for the comparison. The optimization scheme can be expressed in a 
way which is independent of the type of fiuctuation spectrum and of the 
size of the simulations. 



1. Lagrangian Perturbation Theory and Optimization 

For an in-depth discussion of the Lagrangian perturbation approach relevant to 
this work see (Buchert 1994), and for the optimization approach (Buchert et al. 
1994, Melott et al. 1995, Weifi et al. 1995) and references therein. 

Denoting comoving Eulerian coordinates by q and Lagrangian coordinates 
by X, the field of trajectories q = F(X , t) up to the second order on an Einstein- 
de Sitter background reads: 

F = X + q^{a)Vx^^^\X) + q2{a) V x^^^\X), (1) 

where the time-dependent coefficients can be expressed in terms of the expansion 
function a{t) = {j-fl^: 
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The perturbation potentials have to be constructed by solving iteratively the 
two boundary value problems: 

Ax*(') = I{^,^,k)t^. (4) 

Ax^'^ = 2II{^>% (5) 

where / and // denote the first and second principal scalar invariants of the 
tensor gradient {^^]\)- Under the restriction of parallelism of the initial peculiar- 
velocity field and the peculiar-acceleration field, we can set $(^) = being 
the initial peculiar- velocity potential), and the first-order fiow-field (1) reduces 
to Zel'dovich's approximation, as discussed in (Buchert 1994) and references 
therein. 

The optimization scheme used here for the Lagrangian perturbation theory 
was discussed in Coles et al. (1993) and further investigated for higher orders 
in Melott et al. (1995). The idea is to smooth away some of the small-scale 
fiuctuations in the initial data. In a previous work, Melott et al. (1995), we found 
that convolution of the initial density field with a Gaussian window function 
W{k,kgs) = exp( — k^/kg^) of appropriate width kg^ considerably improves on 
the performance of the Lagrangian perturbation schemes. Specifically, here we 
run a search-loop over different smoothing lengths applied to the initial data of 
the first and second order Lagrangian approximations. Then we compare the 
resulting density fields with those obtained with a PM-code (Kates et al. 1991) 
by determining the cross-correlation coefficient S = {61^2) I <^i<^2 of the two fields 
(here, the mean (. . .) is evaluated over all cells of the discrete density fields). As 
optimal we define that smoothing length kg^ which maximizes S . 

In a further analysis, we compute a scale-dependent cross-correlation S{Rg) 
by smoothing the two final density fields ^1, 62 with various filter widths Rg. 

2. Numerical simulations with COBE-normalized SCDM and BSI 
initial data 

In order to follow the nonlinear evolution of the formation of structure we have 
performed N-body simulations using a standard PM code (Kates et al. 1991) 
with 128"^ particles on a 256"^ grid (Kates et al. 1995). The universe is assumed 
to be spatially fiat (i7 = 1). It is dominated by cold dark matter. We consider 
here simulations which were performed in boxes of 500 Mpc, 200 Mpc 
and 75 Mpc. The simulations were started with the power spectrum P{k) 
of density perturbations calculated si z = 2b, 

P{k)={^-l)\^ = \{^f^\k)T\k), (6) 

where = 2H~^ denotes the horizon scale. The primordial perturbation 
spectrum $ is either the Harrison- Zel'dovich spectrum ($ = const) of the SCDM 
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model or the spectrum with broken scale invariance calculated from a double 
inflationary model (Gottlober et al. 1991). In the BSI model the primordial 
spectrum is of Harrison- Zel'dovich type both in the limit of very large and very 
small scales. 

We have used the CDM transfer function T(k) of Bond and Efstathiou 
(1984). For all simulations we have normalized our spectra using the 10°-variance 
of the CMB fluctuations ax = (30±7.5)//K of the flrst year COBE data (Smoot 
et al. 1992). Thus, the power spectrum of the BSI model shows less power on 
small scales than the SCDM model. The scales of nonlinearity deflned in eq.(7) 
below are Afp^^ = 27/1"^ Mpc and X^f^ = Ih'^ Mpc at the time ^ = 0. In 
Fig. 1 we show the linear BSI and SCDM spectra and indicate the box sizes and 
the resolution of our simulations. 




Figure 1. The linear power spectra of 



0001 k [h/;,;,^] our SCDM (dashed line) and BSI (solid 

line) models at time z = Q. 

The flrst and second year COBE data were analyzed by many authors using 

different statistical techniques. The new normalization of Gorski et al. (1994) 

is about 25 % higher than the normalization which we used in the simulations. 

Consequently, in this normalization the scales of nonlinearity would increase to 
\SGDM ^ 33/^-1 ^BSI ^ ^p^_ 



3. Results and Discussion 

In an analysis of both models at different times (z = 0,1,2) we flnd that the 
optimal kgs can consistently be estimated from the scale of nonlinearity kni of 
the considered spectrum and time evolution, kni being deflned by 

As can be seen from Fig. 2, kg^ — lAbkni for the flrst-order and kg^ — l-lk^i for 
the second-order Lagrangian perturbation solutions. It can be seen from the lines 
in Fig. 2, which are regression flts to the scatter plot of kg^/k^i versus n(kni), 
that the optimal smoothing length does practically not depend on the slope 
n(kni) of the power spectrum at the nonlinearity scale, and is thus independent 
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Figure 2. The optimal 
smoothing length kg^ in 
units of kni versus the slope 
n(kni) of the power spec- 
trum at kni for the BSI and 
SCDM models. 

of the exact shape of the power spectrum. While the smoothing length for 
first-order shows some random scatter with n(kni), the optimal smoothing of 
the second order is much more robust against variations in the initial power 
spectrum. 

The scale-dependent cross-correlation function S(Rg) measures primarily 
whether mass is moved to the right place. Fig. 3 shows its value for the com- 
parison of the density fields at time z = 0. The performance of the optimized 
Lagrangian perturbation schemes is very good down to scales of 1 . . .2/i~^Mpc 
for both the SCDM and BSI models. This quality can clearly not be met by 




Rg(BoxsIze/1 28) Rg(Boxsize/1 28) 

Figure 3. Scale-dependent cross-correlation function S(Rg) of first-order (dash- 
dotted lines) and second-order (dashed lines) optimized Lagrangian perturbation 
theory as well as of "chopped linear theory" (dash-dot-dot-dotted lines) corre- 
lated with the results of the numerical simulations for the SCDM (left) and 
BSI (right) models, evolved to 2; = 0. Thin lines are for the 500/i~^Mpc boxes, 
medium thick lines for 200/i~^Mpc and thick lines for the 75/i~^Mpc boxes, 
respectively. 
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the Eulerian linear theory. (For this comparison we use an improvement on the 
Eulerian linear theory as proposed in Coles et al. (1993), called "chopped linear 
theory"; it complies to both Sdin > —1 and {Sdin) = 0, even at late times. In 
terms of the linearly evolved we set 1 + 6 dm = a{l-\- Sun) if ^Un > — 1 and 
otherwise, where a is a normalization constant keeping the total mass the same.) 

The two-point correlation function ,^(r) of the numerical simulations is re- 
produced well by the optimized Lagrangian perturbation schemes down to about 
the "correlation length" ro, ■^(j'o) = 1 (Fig. 4). Below ro, which itself is under- 
estimated by about 10%, the amplitude of ,^(r) of the optimized Lagrangian 
perturbation schemes drops well below that of the numerical simulations. Still, 
the use of our optimization scheme leads to a considerable improvement in ,^(r) 
on small scales, especially for models with a large amplitude of small-scale fluc- 
tuations, like the SCDM model. 



SCDM z = BSI z = 




r [h-' Mpc] r [h- Mpc] 

Figure 4. Two-point correlation function for the SDCM (left) and BSI (right) 
models at 2; = 0. Thick solid lines are for the numerical simulations, the dotted 
line is the prediction of the Eulerian linear theory, the thick (thin) dash-dotted 
and dashed lines are for first- and second-order optimized (unoptimized) La- 
grangian perturbation solutions, respectively. 

The optimized Lagrangian perturbation schemes reproduce the density field 
very good down to the scales of groups and clusters of galaxies, even deep into 
the nonlinear regime (i.e. today). Their inherent advantage is their high speed 
of execution, their complexity is comparable to one time step of a conventional 
particle-mesh N-body code. Thus, they allow the quick production of a large 
ensemble of simulations, e.g. for the evaluation of different cosmological models 
or for doing statistics over an ensemble of different realizations of a single model. 
Furthermore, the analyticity of the mapping (1) allows not only the local deter- 
mination of the resulting shear stresses, e.g. for the introduction of local biasing 
schemes, it can also be used to enhance the particle density by interpolation of 
the mapping (1), so that one can use statistics with a high selection effect, like 
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the simulation of galaxy catalogues with geometry and luminosity selection, even 
for very large-scale realizations (on Gpc-scales). The analysis of the singularity 
structure of the mapping (1) in order to identify building blocks of large-scale 
structure offers a way to analytically relate initial data to present-day structure 
(Buchert et al. 1995). 
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